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We develop the Hartree-Fock-Bogoliubov theory at finite temperature for Bose gas trapped 
in the two dimensional optical lattices. The on-site energy is considered low enough that 
the gas presents superfluid properties. We obtain the condensate density as function of 
the temperature neglecting the anomalous density in the thermodynamics equations. The 
condensate fraction provide two critical temperature. Below the temperature Tc\ there is 
one condensate fraction. Above two possible fractions merger up to the critical temperature 
Tc 2 - Then the gas provides an first order transition at temperature above Tq 2 where the 
condensate fraction is null. We resume by a finite-temperature phase diagram where can be 
identify three domains: the normal fluid, the superfluid and the superfluid with two possible 
condensate fractions. 
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I. INTRODUCTION 


The Bose-Einstein condensate trapped in the periodic potential provided a peculiar phe¬ 
nomenon! in which there exist the quantum phase transition in the ultra-cold atoms similar 

n q 

to the phase transition seen in crystals 1, [2j]. In this transition, the variation of external pa¬ 
rameter related to the trap potential modifies the non-local superfluid properties into the local 
properties of the Mott insulator phase. 

At low temperatures, the superfluid-to-Mott-insulator (SF-MI) quantum phase transition was 
proposed in 1998 by Jaksch et al. jii] for ultra-cold gases. They considered a Bose-Hubbard model 
to describe the Bose gas on the optical lattices, and they showed the quantum phase transition 
by changing the on-site energy, the hopping matrix element or the chemical potential. Four years 
later, Greiner et al. 4j] detected the SF-MI phase transition by imaging atoms released from optical 
lattices with different periodic potentials. Above a critical on-site energy, the atomic distribution 
presents the correlation function concerned to the Mott-insulator. In this phase, the position 
of atoms become a tool to probe and measure thermodynamic properties such as pressure, 
temperature, and transport properties jf]]. 

Another way to observe the SF-MI phase transition is analysing the collective excitations in 
the gas trapped in optical lattices. The low-lying excitation energy of the gas is achieved in ID 
via lattice modulation [7] or by Bragg spectroscopy {s] . In the superfluid (SF) phase, the energy 
spectrum presents a gapless behaviour which can be explained within a Bogoliubov theory jij, 


or 


in terms of path integral formalism [10] or by Green’s function formalism 


111 - 


In Mott insulator (MI) phase, the gapped excitation mode is found by the mean-field theory 
9] using the particle-hole excitations frame of the superfluity modes. The mean-field theory in 
9| neglected the quantum fluctuations, then the description of the system near to SF-MI phase 


transition is unavailable and, therefore, not confirmed the gapped mode 


12], The SF-MI transition 


theory which treats the kinetic energy as a 


is analysed in terms of a strong-coupling perturbation 
mean-field and deals exactly the on-site repulsion [3, 13j|. This approximation gave a estimative 
for the SF-MI transition limits. 

The thermal gapless properties of ultra-cold bosons in optical lattices were discussed by Gerbier 
in the ref. [14]. However, his analysis were done considering only the Mott insulator phase. At 
superfluid phase, a finite-temperature study of ultra-cold bosons in 2D optical lattices were done 
in [3] by quantum Monte-Carlo simulation taking account the effects of harmonic confinement. 
In this reference, they observed the spatial inhomogeneity allows the coexistence of three phases: 
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the usual ground state superfluid, the Mott insulator phase and the normal fluid phase, 
simulations can be compared with the phase diagram provided by the experimental results 


t: 


reir 



The phase diagram is describe as a function of the temperature and the interaction strength 
whose definition is the ratio between the on-site interaction strength and the hopping matrix 
element. The phase diagram for the homogeneous 3D system at unity filling is plotted using the 
quantum Monte Carlo simulation in the reference [17] and shown experimentally in [18]. At zero 
temperature, the system undergoes a SF-MI quantum phase transition at the critical interaction 
strength. At temperatures different from zero and near to the critical interaction strength, they 
observed a suppression of the critical temperature where the condensate fraction is null and all 
particles are in the thermal cloud. 


In this work, we investigate the thermal behaviour of the ultra-cold Bose gas in 2D optical 
lattices on the superfluid phase and far from the SF-MI phase transition by a finite-temperature 
mean-field theory. This theory is known as Hartree-Fock-Bogoliubov (HFB) theory and departing 
from the Bose-Hubbard model. Although the HFB theory provide a gapped spectrum even far 
from the quantum phase transition, we avoid this problem by the Popov approximation where the 


anomalous densities are neglected ad hoc , 


m 


and experimentally observed in [l2], the condensate fraction has the superfluid properties 


dence the spectrum is gapless and, as theoretically shown 


and depends on the temperature and the interaction strength. Then, in HFB theory, the superfluid 
phase has, as order parameter, the fraction of condensate. We determine the order parameter as 
function of the temperature at a regime below to the quantum phase transition. 


In sec. QH we obtain the HFB theory considering the superfluid state as a vacuum of the 
quasi-particles at finite temperature, we determine the limits to the Bogoliubov theory and we 
implement the Popov approximation to derive the condensate fraction. The numerical solution for 
the condensate fraction in the Popov approximation is evaluated in the sec. IIII1 Finally, in sec. IIVI 
we discuss the properties of the system near the thermal phase transition according to the solution 
of the mean-field theory. 
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II. HARTREE-FOCK-BOGOLIUBOV APPROXIMATION AT FINITE TEMPERATURE 

The dynamics of the 2D ultra-cold boson gas in the optical lattices can be done by the second 
quantized Hamiltonian, 

H = J ^(x) + ^o(x) + V T (x)) V’(x)d 3 x + (1) 

+ 1 ^ a sh f ( x )-0t ( x )-0( x )-0( x )(i 3 X 

2 m J 

where V’(x) is the boson field operator corresponding to the atom with the mass equal to m, 
Vt(x) = f Ej = i is the magneto-optical trap and Vo( x ) = Vo(sin 2 (7nc/a) + sin 2 (7n//a)) is 
the external periodic potential formed by the intersecting laser beams with Vo being the trap 
depth proportional to the laser intensity, a being the less distance between two sites. In this 
Hamiltonian, the inter-atomic interaction corresponds to the low-energy collision characterized by 
the .s-wave scattering length, a s . 

We expand the field operators in Wannier functions whose asymptotic behaviour is adjusted 
by the Gaussian function, w(x) = e“( x2+J/2 )/ 2 <t2 /(7t 1 / 2 <7) 1 / 2 with a = h/ {rYVoEr) 1 ^ and Er = 
h 2 ir 2 /(2ma 2 ) is the recoil energy. Considering only the lowest vibrational state in each site, the 
field operator can be written as, V’(x) = Yh w(x — Xj)aj, where x* is the position of the site labelled 
by i in the 2D lattices with N s sites. Thus, in terms of the Wannier operators a*, the Hamiltonian 
become a Bose-Hubbard (B-H) Hamiltonian, 

H = -tJ2 a i a j + \ U “I' a i a i a i’ ( 2 ) 

(i,j) i 

where t = — / w*(x — x t ) + Vo( x )^ w(x — Xj)d 3 x is the hopping parameter and U = 

f |rc(x)| 4 ci rf x is the on-site interaction strength. The hopping parameter is the same for any 
site due to neglect the magneto-optical trap effects, Vr(x) = 0. 

In order to identify the condensate state, we expand the Wannier operators in the 2D quasi¬ 
momentum space, a t = -^^^'^2] i a\ i e~ l ^ kxXi+kvyi \ where the component of quasi-momentum, k, 
assume discrete values give by kj = 2irri/(Mja) with = 0, ±1, ±2, • • • , ±(Mj — l)/2 and Mj is 
the even number corresponding to the number of site in the j-direction. 

The condensate state is a coherent state of No single particles with quasi-momentum equal 
to zero. We represent this state in the B-H Hamiltonian by shifting the zero quasi-momentum 



5 


operator, Uk = ^k + ^/NqSq k- Then we have the shifted B-H Hamiltonian as, 


H — (—eo — h)Nq + + v // ^o(—eo — ^ + gNo){bo + 6q) + 


( 3 ) 


+ ^(~ e k — P + 2<7-/Vo)& k &k + 2^o(^k&—k + ^k^-k) 


+ 


+ z 0 g^2(bl +(i blb q + 6j c _ q 6 k ^ q ) + f ^ &k+q & l- q bkbk '> 

kq kk'q 

where ek = 2t(cos(k x a ) + co s(k y a)) is the energy of the free quasi-particles, g = U/N s is the 
interaction parameter and p is the chemical potential in which ensures a well-defined particle 
number even for the particle reservoir contact. 

For the purpose of decoupling normal modes, we introduce the Bogoliubov transformation, 
6 k = '“kCk + HcC-k where and «k are the real Bogoliubov coefficients and obey u k — u k = 1 to 
keep the commutation relations. The Bogoliubov coefficients are adjusted for becoming the B-H 
Hamiltonian in diagonal form, H = E$ + £fcC k Ck, where Eq is the ground state energy and Ek 
is the excitation energy. 

Assuming the time invariance of the Bogoliubov coefficients, we derive the excitation energy 
as, Ek = \Jh\ — A 2 , and the Bogoliubov coefficients are, v\ = \ — lj and UkVk = — 

where the parameters h k = e o — Q; + <?(Ao — Y2k K k) and A = g(No + Yhk K k) depend on the 
variational parameters, pk = v{k,T) + (1 + 2v(k,T))v%. and Kk = (1 + 2n(k,T))u k Vk^ These 


parameters are thermal averages of shifted operator and depend on the temperature 19]. The 
function v(k,T) = — l) -1 is the Bose-Einstein distribution of a quantum particle with 

energy equal to Ek■ 

Finally, we can obtain the chemical potential in terms of the number of the condensed particles, 
Nq, the depleted particles, YlkPk and the paring, Ylk K k- The chemical potential is written as 
P = — + g[N 0 + + K k )]• 


A. Particle number 


To evaluate the condensate fraction no = Ao/lV s as an explicit function of the dimensionless 
interaction strength, U/t, and the temperature, T, we introduce the total density in each site by, 
n = ]Ck( a k a k)/-^s- I n terms of the Bogoliubov parameters, the depletion density is given by, 
n — no = X^k Pk/^Ys- Hence, we have the condensate density in each site as, 


1 

nn = n - 

N„ 





( 4 ) 
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where, /ik = Cfc + U{n o — mo), A = U{uq + mo) and Cfc = 4i[sin 2 (fc x a/2) + sin 2 (fc y a/2)]. 

The anomalous density in each site, mo = CCfe K k)/N s , can be an explicit form that depend 
on the condensate density, mo = — n oCCfc co ^ a -(Ek/2k,BT)/E^)/{2,/U + cot\\(Ek/2kBT) / E^). If 
the anomalous density is different from zero, we have a gap in the excitation spectrum given by, 
lim fc _x) E k = 2U (—n 0 m 0 ) 1//2 . 


B. Popov approximation 


In the HFB theory, the excitation spectrum presents a 


therefore, not satisfies the Hugenholtz-Pines theorem 


gap 


20 


for long-wavelength excitations and, 


21fl. Indeed, the gap in the spectrum 


for long-wavelength can not explain the superfluid properties, as the absence of viscosity [22. |23(]. 
widely identified due to the phonon behaviour of the excitation spectrum in ultra-cold alkaline 
atoms 24]. 


At very low temperatures and in the very weakly interaction regime, the condensate fraction 
is high enough to neglect the depletion, no ~ n, and the anomalous density mo = 0. By this 


approximation, known as Bogoliubov approximation, the excitation energy, E\ = ^/Ck(Ck + 2 nil), 
is gapless for the long-wavelength limit and the condensate density depends only on the total 
density, n. 

Otherwise, we can neglect ad hoc only the anomalous density, mo = 0. The so-called Popov 
approximation provides a non-linear integral equation to the condensate density given by, 


riQ = n 


1 

w. 


E 


coth (E k /2k B T) (Cfc + n 0 U) - E k 
Ek 


(5) 


with a gapless spectrum equal to E). = \J (k(.Ck + 2rioU). 


III. NUMERICAL RESULTS 

We obtain the condensate fraction, no/n, as function of the temperature and the dimensionless 
interaction strength parameter, U/t, solving numerically the non-linear equation ([S]) derived by the 
Popov approximation. 

At zero temperature, we plot in fig. |T] the condensate fraction as function of the interaction 
strength for four sizes of the 2D squared lattice. When we increase the size of the lattice, the 
condensate fraction converge to an asymptotic curve regarded to the thermodynamic limit. Then, 
as the difference between the curves provide for lattice with 101 x 101 and 501 x 501 is small and 
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for the sake of the computation time, we assume the size of lattice equal to 101 x 101 can be 
considered in the thermodynamic limit. 

In the ref. 9], they analysed the limits for the strongly interacting regime, U/t —>• oo, at zero 
temperature and they concluded that the quantum phase transition could not be identified in the 
Bogoliubov approximation. However the arguments used to justify the absence of phase transition 
are based on the Popov approximation indeed. We observe the same result in the fig. |Tj for any 
size of lattice using the Popov approximation. Differently from them, it is not clear for us the rule 
of the Bogoliubov approximation in this lack of the quantum phase transition. 



FIG. 1: The condensate fraction no/n in two-dimensional lattice as function of the dimensionless interaction 
strength parameter UIt at zero temperature. The graph exhibits the approach to the thermodynamic limit in which 
the ration N/N s = n remains constant. This approach is realized when the size of the lattice exceeds 10 4 sites. 


In order to obtain the temperature-dependent condensate fraction, we introduce the function, 


f( x ) = i-x- — -y- 

v ’ nN q ^ ? 


1 v^l ( coth ( \/tk(ek + 2 nUx)\ e k + nUx 


\ J y^kCek + 2 nUx) 


- 1 


( 6 ) 


where the variable x running from 0 to 1. The function f(x) is determined by the temperature 
T and the interaction strength parameter U/t and its roots are the condensate fractions at this 

















X 


FIG. 2: The function f(x) (eq. 0 is plotted for different values of the temperatures and for n = 1.0, U/t = 10.0 
and Ns = 101 x 101. This graph illustrates the numerical method used to find the roots of ,f(x). 


temperature and for this interaction strength. In fig. [2] we plot this function for U/t = 10.0 and 
at different temperatures. The plot exhibits three different types of the solution for equation ([5]) . 
There is one solution for 0 < T < Tci, where Tci is the temperature in which the function f(x ) 
presents a root at x = 0. When we increase the temperature, we find two solution which merger 
at T = Tc 2 - At temperatures larger than Tq 2 , there is no solution and, therefore, there is no 
condensate fraction in the lattice. 

In the fig. [3] we plot the condensate fraction as function of the temperature for different values 
of the interaction strength. In the non-interacting case U/t = 0.0, the fraction decrease when the 
temperature increase and reach a null fraction, i. e., the condensate disappear at the temperature 
Tci = 1.6954 (in units of ks/t). Above this critical temperature the gas is entirely a thermal Bose 
fluid and the quasi-particles have free energy dispersion in the limit of the large wavelength. 

For the interacting case U/t ^ 0.0, there are only one condensate fraction in the interval of the 
temperature 0 < T < Tc i- In this condition, the values of the fraction decrease as temperature 
increase. For the temperature in the range of Tci < T < Tc 2 , there are two possible condensate 









FIG. 3: The condensate fraction no/n as a function of the temperature T for different values of the interaction 
strength parameter U/t in the same conditions of [5] Three interacting regimes are deal in the figures. The 
non-interacting gas is describe by the solid curve, the weakly interacting gas is shown in the traced line, and the 
strongly interacting gas is represented by the dotted curve. 


fractions. Differently from the largest fraction which decreases as the temperature increases, the 
smallest fraction increases up to merger the highest fraction at the second critical temperature Tc 2 - 
At temperature above 7c2, there is no solution for f(x) and, therefore, there is no condensate in 
the lattice. 

The fig. [3] shows that exist an interval of the temperature in which the strongly interacting 
Bose gas (U/t = 10.0) presents the condensate fraction larger than the condensate fraction of the 
weakly interacting Bose gas (U/t = 0.1). This unusual behaviour is contrary to what intuition or 
common sense would indicate. 

The phase diagram shown in the fig. [4] has three regions correspond to the possible values of the 
condensate fraction. The region SF where there is one condensate fraction, CSF in which coexists 
two possible condensate density and the region NF composed by the thermal free quasi-particles. 
As we increase the interaction strength parameter, the CSF region increases monotonically. 










FIG. 4: Finite-temperature phase diagram as a function of dimensionless parameter U/t and temperature T of 
the Bose gas trapped in two dimensional lattice with size of Ns = 101 x 101 and density equal to n = 1.0. 


Finally, in the fig. [5] we present the condensate fraction no /n in two dimensional lattice as 
a function of the interaction strength parameter at several temperatures. At temperature below 
the critical temperature Tqi, the condensate fraction increase when the interaction parameter 
increase. This growth is provided by the positive contribution of the interaction. The repulsive 
interaction promotes the non-localization of the atoms in the lattice. Although the effect of the 
interaction, the entropy plays the leading role for large values of the interaction parameter. Then, 
the maximum value of the condensate fraction is reached around U/t ~ 0.5, and the fraction 
decreases monotonically up to zero in the limit of strongly interaction. At temperatures higher 
than the critical temperature Tq 2 , the curves present a minimum value to the interaction parameter 
in which the two possible solution of the condensate fraction rise. As the case of low-temperature, 
there is a balance between the interaction and the entropy. For interaction parameter close to the 
critical parameter, the role of the interaction is more determinant than the entropy. However the 
entropy dominates for the limit of the strongly interacting regime. 





FIG. 5: The condensate fraction no/n as a function of the dimensionless parameter U/t at different temperatures 
(in units of t/fes) with n = 1.0 and Ns = 101 x 101. At temperatures below Tbi, the curves are positive functions 
with only one maximum point. At temperature above, the curves present a minimum critical value for the 
dimensionless parameter Uft where below there is no condensate at this temperature. For above the critical 
dimensionless parameter, the condensate can coexist in two possible condensate fraction. 


IV. SUMMARY 

We develop the Hartree-Fock-Bogoliubov theory at finite temperature for Bose gas trapped in 
the two dimensional optical lattices. This mean-field theory presents anomalous pairing density in 
which provides a gap in the energy spectra for large wavelength excitation. 

To obtain a gapless excitation energy, we consider the Popov approximation where the anoma¬ 
lous density is neglected ad hoc. In the Popov approximation, we derive a non-linear equation 
whose solutions give us the condensate fraction in the system. Considering the size of the lattices 
large enough to be in the thermodynamic limit, we obtain curve of condensate fraction as function 
of the temperature. These curves present two critical temperature, one of them Tq\ define an 
interval where we have a single order parameter. Above this temperature, there are two conden¬ 
sate fraction that merger at the critical temperature Tq 2 - Then the gas provides an first order 
transition when temperature pass over Tq 2 and, in this condition, the condensate fraction is null. 












12 


We identify an interval of the temperature in which the strongly interacting gas has the condensate 
fraction large than the weakly interacting gas. This unusual effect is due to the balance between 
the role of the entropy and the interaction. 

A phase diagram as a function of interaction strength U/t and temperature T is done calling 
attention to the regions where occur the normal fluid phase, the superfluid phase and the superfluid 
phase in which coexist two different condensate fractions. There is a decay behaviour when the 
interaction strength approach to the zero as observed by [17] using the quantum Monte Carlo 
simulation for 3D lattice. Then the HFB theory is qualitatively usefull to describe the condensate 
fraction in the weakly interacting regime. The quatitative comparison between the HFB theory and 
the Monte-Carlo simulation is out of scope and it will be done in future work. Otherwise, for large 
interaction strength, the thermal phase transition temperature increase and, therefore, frustrates 
any possibility to identify the quantum phase transition by the mean-field theory. Indeed, the trial 
function to develop the HFB theory is not appropriated to describe the strongly correlated lattice. 
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